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ABSTRACT 

We present the discovery and a detailed multi-wavelength study of a strongly-lensed luminous 
infrared galaxy at z=0.816. Unlike most known lensed galaxies discovered at optical or near-infrared 
wavelengths, this lensed source is red, which the data presented here demonstrate is due to ongoing 
dusty star formation. The overall lensing magnification (a factor of 17) facilitates observations from 
the blue optical through to 500 /rm, fully capturing both the stellar photospheric emission as well as the 
re-processed thermal dust emission. We also present optical and near-IR spectroscopy. These extensive 
data show that this lensed galaxy is in many ways typical of IR-detected sources at z ~ 1, with both 
a total luminosity and size in accordance with other (albeit much less detailed) measurements in 
samples of galaxies observed in deep fields with the Spitzer telescope. Its far-infraed spectral energy 
distribution is well -fit by local templates that are an order of magnitude less luminous than the lensed 
galaxy; local templates of comparable luminosity are too hot to fit. Its size (D ~ 7 kpc) is much larger 
than local luminous infrared galaxies, but in line with sizes observed for such galaxies at z ~ 1. The 
star formation appears uniform across this spatial scale. Thus, this lensed galaxy, which appears 
representative of vigorously star-forming z ~ 1 galaxies, is forming stars in a fundamentally different 
mode than is seen at z ~ 0. 

Subject headings: galaxies: high-redshift — galaxies: evolution — gravitational lensing 


1. INTRODUCTION 

Strongly gravitationally lensed galaxies offer a unique 
view of the distant Universe, by providing a surfeit of 
photons and access to physical scales not resolvable when 
studying galaxies that have not been lensed. Further- 
more, lensed galaxies are likely to probe a mo re typi- 
cal portion of the luminosity function (e.g., Bayliss et al. 
2010) than strictly magnitude-limited blank held studies. 
The typical bright lensed source, from samples selected at 
optical wavelengths, is a vigorously sta r-forming galaxy 
at z~2 (Bayliss et al. pOTT ; Bayliss Ml with blue colors 
and modest reddening fWuyts ct al.|2012p. The brightest 
of these lensed sources, w hich are typically magnified by 
factors of 20 or more (e.g. , Wil liams fe Lewis|1996 Allarn 
et al.|2007 Sharon et al.fjOl^ — are of particular interest 
in part because they enable studies in wavelength regimes 
where current facilities lack the sensitivity to reach high 
redshift galaxies (see the discussion in Siana et al. 2009). 
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The brightest lensed sources reported to date from dis- 


prototype MS1512-cB58 (Yee et al. 1996 

I to more recent 

discoveries (e. 
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Belokurov et al. 2007 

Lin et al. 2009 

Koester et al. 
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L0D. The reddest ot these is the Cosmic 

Eye at z = 3.L 

74 (]Smaii et al. 2007 Siana et al. 2UU9). 


limit has pro vided a sensitivity boost for sub-mm g alaxy 
surveys (e.g., Smail et al. 1997 Egami et al. 2010), but 
typically these sources are magnified only by mns of per- 
cent, with no multiple imaging. As such, there are only 
a small number of robustly confirmed gravitati onally- 


lcnsed sub-mm galaxies with large magnifications (Rigby 
et al.||2008t ISwinbank et aLl|20101 ILestrada et al. 1 1 1 


Frayer et al.||201~ Combes et al. |[2012 Fu et all 2 0 1 2 [ ). 
Expectations are that the new generation of millimeter 
and submillimeter telescopes should detect large num- 
bers o f lensed galaxies at t hese wavelengths (e.g., Her- 
schel: |Negrello et al.||20T0} [Gonzalcz-Nuevo et al.||2012| 


The South Pole Telescope: Vieir a et al.||2010p T~ 

In this paper, we report the discovery of an unusually 
red and bright lensed galaxy, at a redshift of z = 0.816. 
This galaxy is multiply imaged and magnified by a to- 
tal factor of 16.8^6 2 by a foreground galaxy cluster at 
z = 0.237. We present photometric measurements of the 
source from the restframe ultraviolet to the far infrared, 
as well as both optical and near-infrared spectroscopy. 
From these data we construct a lens model, which we 
use to invert measurements to unlensed quantities in the 
source plane when required. These data show that the 
source is a luminous infrared galaxy (a LIRG) that is 
typical of star-forming galaxies at that epoch. §2 de- 
scribes how the lensed source was discovered as part of a 
larger effort to find and characterize lensed galaxies. §3 
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describes the follow-up observations, as well as the basic 
measurements of the source. §4 details the lens model 
and properties of the lensed source in the image plane; 
these measurements are discussed and contexualised in 

§5- 

Throughout the paper we adopt a flat cosmology with 
f 1m = 0.3 and H 0 = 70 kms -1 Mpc” 1 . All magnitudes 
are reported on the AB system unless otherwise noted. 
All on-sky images show North up and East to the right. 


2. DISCOVERY OF A RED HERRING AND A RED LENSED 

GALAXY 

The cluster SDSSJ 1438+1454 was selected as a pos- 
sible strong lensing system (M. D. Gladders et al., in 
preparation) from SDSS imaging, likely due to the pres- 
ence of a blue arclet feature 1875 from the brightest clus- 
ter galaxy (BCG), with color, geometry with respect to 
the BCG, and a hint of curvature all suggestive of lens- 
ing, as shown in Figure 1. As described in Appendix A, 
follow-up spectroscopy from the ARC 3.5m showed that 
this initial arc candidate is not lensed, but is rather a 
cluster galaxy at z = 0.231. 

However, another arc-like feature, this one faint and 
redder than the BCG, was just visible in the SDSS imag- 
ing. We had obtained warm Spitzer observations with 
IRAC at 3.6 and 4.5 /tm for a larger sample of ~ 100 lens- 
ing systems through Spitzer program 70154 (PI: M. D, 
Gladders). This Spitzer imaging (Figure 2, and Table 1) 
clearly showed that this source was infrared-bright, very 
red in the g — 3.6 /jm and g — 4.5 /mi colors, and likely to 
be a multiply imaged lensed galaxy. This source is the 
subject of this paper. Fo llowing the naming conventions 
established in Koester et al. (2010) we refer to the lens- 
ing cluster as SDSSJ 1438+1454 and the lensed source as 
SGAS 143845.1+145407 (SGAS 1438 in short hereafter). 
The source name coordinates are derived from the posi- 
tion of the most isolated complete image, which is 5.6" 
almost due south of the BCG. 


3. FURTHER DATA AND MEASUREMENTS 

3.1. Optical and Near-IR Imaging 

Motivated by the identification of 
SGAS 143845.1+145407 as a possible lensed source 
from the SDSS and Spitzer imaging, we observed the 
field with other facilities. Near-IR imaging in the YJH 
and K s bands was acquired during commissioning of the 
FourStar camera at the Magellan Baade 6.5m in early 
2011. These near-IR images robustly confirmed that 
this was indeed a lensing geometry, as shown in Figure 
3, and further demonstrated that the lensed source is 
extremely red. The FourStar JH and K s images were 
flux-calibrated directly to th e Two Micron All Sky 
Survey point source catalog (|Skrutskie et al.| [2006 ), 
which is possible because the field of view of FourStar 
is large; calibration uncertainties are a few percent in 
these bands. The Y-band image was calibrated using a 
standard star observatioif^l acquired immediately prior 
to the science observation at a similar airmass. This 
calibration agrees with repeated Y-band calibrations 


10 Calibrations of the Y-band and other medium-band filters 
available in FourStar h ave been undertaken in support of the Z- 
FOURGE project |Spitler et al. | (|2011| and will be published in 
detail elsewhere. 


of the instrument over months of observing to within 
0.015 magnitudes; we thus consider that the Y-band 
data have calibration uncertainties of a few percent. In 
all four bands these uncertainties are sub-dominant to 
other sources of error. 

Optical imaging in griz was obtained with the GMOS 
imaging spectrograph on the Gemini North 8.2m tele- 
cope; zeropoints and color terms were established by 
comparison to SDSS-measured stars in the same field and 
are uncertain by a few hundredths of a magnitude. 

The GMOS/Gemini-N and FourStar/Magellan im- 
agery is shown in Figure 3, with images centered on the 
brightest cluster galaxy. Table 1 summarizes the details 
of the observations and data. Photometry of the lensed 
source is reported in Table 1, with details of the measure- 
ment process including treatment of the central galaxy 
given in §3.4. 


3.2. WISE photometry 

SGAS 1438 was included in the Wide-Field Infrared 


Survey E xplorer (WISE) preliminary release (Wright et 
aLpOlO). As it is extended in the WISE bands, we per- 
form "aperture photometry rather than take the profile- 
fit photometry from the WISE preliminary release source 
catalog. Though the source is detected in the bluest two 
bands of WISE, these data are completely superseded by 
the higher spatial resolution and deeper imaging from 
Spitzer-IRAC, so we report here only the WISE mea- 
surements at 12 /tm and 22 /tm. These wavelengths are 
uniquely covered by WISE for this source, and are impor- 
tant since they measure the strength of rnid-IR aromatic 
emission features. The details of this data and measured 
photometry are reported in Table 1; details of correc- 
tions for contamination by light from the BCG are given 
in §3.4. 


3.3. Herschel photometry 

We obtained mid- and far-IR photometry for 
SGAS 1438 in discretionary time on Herschel using the 
PACS and SPIRE instruments, with a total clock time 
of 2.1 hr. Details of the observations are given in Table 

1 . 

For PACS, we performed aperture photometry on the 
level 2 pipleline-produced science images, using annular 
apertures with radii of 10,12,16 pixels at 70 /tm, 12,14,16 
pixels at 100 /tm, and 10,12,16 pixels at 160 /tm, respec- 
tively. We applied aperture corrections as given in tech- 
nical memo PICC-ME-TN-037 vl.0. We applied color 
corrections as given in technical memo PICC-ME-TN- 
038 vl.0, by interpolating to a 27 K blackbody from the 
blackbodies tabulated at 20 and 30 K (3 These color cor- 
rections were 1.09, 1.00, and 0.97 at 70, 100, and 160 /tm. 
At 70 /tm where the spectral slope is particularly steep, 
the uncertainty on the color correction is ~ 10%, and 
thus contributes significantly to the photometric uncer- 
tainty. 

For SPIRE, we performed aperture photometry on the 
level 2 science images, using annular apertures, following 


11 The PACS photometry is quite sensitive to the choice of color 
correction; we determined by fitting the photometry that a T = 
27 K blackbody was an adequate fit to the raw (not de-redshifted) 
photometry. We experimented with modified blackbodies, where 
prnodified _ ,+ (T) . but found these did not provide a better 

fit. 
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the procedure of the SPIRE Photometry Cookbook. In 
the Rayleigh- Jeans limit of F v oc ^ 2 , the color correc- 
tions tabulated in the Cookbook are 0.94, 0.95, 0.94 at 
250, 350, and 500 /jm. Resulting PACS and SPIRE AB 
magnitudes are listed in Table 1, converted from mea- 
sured flux densities f v to vo.ab = 8.9 - 2.51og (/„), with 
ft, in Jansk ys following the definition of AB magnitudes 
Oke] 1974). 

SGAS 1438 is sufficiently bright in PACS and SPIRE 
that the photometric uncertainty is dominated not by 
photon noise, but by a number of systematic effects: un- 
certainty in the color corection, aperture correction, and 
centroiding, and the fact that the source is somewhat re- 
solved in the shortest-wavelength bands. We estimate 
these effects as contributing a total of 13% uncertainty 
on the flux density in each band. 

At its band with highest spatial resolution, 70 /jm, Her- 
schel partially resolves SGAS 1438, in that it splits the 
southern image from the blended other images. Fitting 
with two GALFIT components (Peng et al. (2010); see 


3.4), we find that the southern image contains 22 ±5% 
the 70 /irn light. 

3 . 4 . Subtracting the Central Galaxy and Final 
Photometry 

The field of SDSSJ 1438+1454 is complex: the mul- 
tiple images of the lensed source encircle the brightest 
cluster galaxy, which contributes significantly to the flux 
in the optical, near-infrared, and Spitzer 3.6 and 4.5 /im 
bands, and possibly the WISE 12 /mi and 22 /im bands 
as well. This complexity requires an atypically involved 
analysis to isolate the flux contribution of the source at 
each wavelength. To do this, we fit the BCG in each of 
the optic al and near-IR images using GALFIT (Peng et 
al. 2010). The appropriate point spread function (PSE) 
model for each image was constructed from fitting point 
sources in that image, with a set of PSF objects fit using 
a multi-component moffat+gaussian model in GALFIT, 
and then scaled and summed at a fiducial center to pro- 
duce the final PSF model for each image. For the BCG 
fits, it was found that a multi-component Sersic model 
- with a small central component that is rounder and 
about 3-4 magnitudes fainter than the dominant com- 
ponent - was needed to achieve an acceptable fit. That 
the BCG requires more than one component for an ac- 
ceptable fit is neither uncommon nor unexpected, and 
we have confirmed that this result is not due to a mis- 
construction of the PSF reference. The lensed source 
was masked iteratively, with components revealed after 
an initial BCG subtraction remasked and the fit iterated 
across all bands with the same mask. Figure 4 shows 
the final BCG-subtracted images for a few representa- 
tive wavelengths. 

The resulting BCG models do not show any signifi- 
cant evidence for color gradients in the BCG, and so 
to subtract the BCG flux from the Spitzer 3.6 /tm and 
4.5 /t m images we took the best fit K s -band BCG model, 
convolved it with a PSF constructed as above for the 
IRAC images, and then subtracted a scaled version of 
the resulting BCG image by eye. We did not use GAL- 
FIT to perform the subtraction automatically because 
the Spitzer-IRAC PSF badly blends the lens with some 
of the images of the source and thus the minimization 
calculation in GALFIT performs poorly absent more ex- 


tensive masking of the lensed source. 

For the IRAC data, and all images at shorter wave- 
lengths, we then measured the magnitude of the isolated 
southern i mage of the source using the techniques de- 
scribed in |Wuyts et al. (2010). We also measured the 
ratio of the total flux of the source to this southern 
image in the four reddest bands, where the flux ratio 
of SGAS 1438 to the BCG is highest. This ratio is 
3.49+0.17. Final magnitudes reported in Table 1 are con- 
structed from the measured magnitudes of the isolated 
southern image rescaled to the total image using this ra- 
tio. The simpler approach of reporting total source mag- 
nitudes in the BCG-subtracted images was avoided due 
to the possible influence of systematics from the BCG 
model subtraction. At bluer wavelengths the BCG is 4-5 
magnitudes brighter than the source, such that aperture 
measurements of the entire source would be completely 
dominated by the BCG; by scaling from measurements 
of the isolated southern image, we trade these potential 
systematics which are difficult to quantify for a noisier 
measurement and a quantified systematic uncertainty in 
the scaling of less than 5%. Note that we do not include 
this systematic 5% uncertainty in Table 1; it is only rel- 
evant when comparing results at shorter wavelengths to 
data from Herschel and WISE, and since the spectral 
energy distribution (SED) fitting (see §4.4) is done inde- 
pendently for the shorter and longer wavelengths it has 
no effect on our conclusions. 

At the long wavelengths sampled by Herschel, the con- 
tribution of the BCG is expected to be insignificant. In 
the WISE 12 /rm and 22 /im bands the contribution from 
the central lens galaxy may be significant, and the poor 
spatial resolution of WISE does not separate the BCG 
from the lensed source. We thus proceed by using an el- 
liptical galaxy SED model, anchored at the reddest wave- 
lengths at which we have sufficient resolution to mea- 
sure the central galaxy light independently {HK and 3.6, 
4.5 /nn). We take as the fiduc ial SED the S pitzer data on 

We use the 


tire bpr 
(2006). 
rlaxv ch 


BCGs in clusters from Egami et al. 
data on central galaxies from the galaxy clusters Abell 
773, Abell 2219, and Abell 2261, all of which are close to 
SDSSJ 1438+1454 in redshift, and are noted in Egami et 
al. as normal central galaxies with no excess flux at long 
wavelengths. From these three g alaxies, and their K , s 
magnitudes in the 2MASS catalog (Skrutskie et al. 2006), 
we compute colors for the 4.5 /tm, 8 /in r, and 24 /im bands 
referenced to the K s and 3.6 /tm bands. Uncertainties are 
taken as the spread in these values across the three galax- 
ies. We measure the K s and 3.6 /tm AB magnitudes of 
the best fitting BCG model from our data, and then use 
the colors derived from the data of Egami et al., with 
simple interpolation, to predict the 12 /tm and 22 /im 
magnitudes of the central galaxy in SDSSJ 1438+1454. 
Final uncertainties are derived from the uncertainty in 
the color, added in quadrature to the difference between 
the predicted magnitudes in the K s and 3.6 /im bands. 
The result is that at 12 /mi the WISE photometry of the 
lensed source is nominally contaminated at 10%±5% by 
the central galaxy, and at 22 /im at 2.5%±1.5%. In both 
cases the uncertainty on this correction is less than the 
photon noise in the measurements. 

3 . 5 . Optical spectroscopy 
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The BCG was observed spectroscopically using the 
ARC 3.5m telescope and the Double Imaging Spectro- 
graph (DIS) on the night of 8 April 2011, to measure 
its velocity dispersion. The 1"5 slit was used, with the 
B1200 and R1200 gratings used on the blue and red chan- 
nels respectively. The gratings were centered at 4750 A 
and 6200 A, which provided coverage from 3870 A-5130 A 
and 6250 A-7390 A and nominal resolutions of R ~2000 
and ~3000 at the wavelength centers of the blue and 
red channels. Three 1800 s integrations were acquired 
under good conditions. Data were reduced to flux- and 
wavelength-calibrated stacked ID spectra using custom 
IDL scripts incorperating procedures from the XIDIpJ 
software package. The resulting spectra have per pixel 
SNRs of typically ~ 10 and ~ 20 in the blue and red 
channels respectively. 


3.6. Near-infrared spectroscopy 

We obtained a spectrum of SGAS 1438 with the NIR- 
SPEC instrument on the Keck II telescope on the night 
of 19 Mar 2011 UT. The weather was clear, with the see- 
ing variable during the night and especially poor while 
SGAS 1438 was observed. The airmass was secz<1.07. 
The slit position angle was 346 degrees E of N, such that 
spectra for the two brightest images of SGAS 1438 were 
simultaneously captured, as shown in Figure 5. Some 
light from the BCG is included in the slit though it is 
spatially distinct from the light from SGAS 1438. We 
used the NIRSPEC-2 filter, low resolution mode, a grat- 
ing angle of 36.48, the 0.76" x 42" slit, and an ABBA 
nod pattern with a 14" nod. This setup covered the 
wavelength range of 1.11 to 1.29 /im. We integrated for 
600 s at each nod position. For telluric flux calibration, 
the star A0 V star HD 131951 was observed immediately 
after the observations, at an airmass of secz=1.06 and a 
distance of 4° from SGAS 1438. 

We reduced the NIRSPEC spectra using G eorge 
Becker’s pipeline as described in Rigby et al. (2011), with 
one exception: since the emission lines were too faint 
to extract a one-dimensional spectrum from each nod, 
we instead summed the sky-subtracted two-dimensional 
frames taken in Nod A, repeated the summation for Nod 
B, and then separately extracted the spectra from each 
summed nod. For the North image we weighted-average 
combined the extracted spectra for each of the two nods. 
For the South image we used only the spectrum from 
Nod B, as the spectrum from Nod A suffered persistance 
from the trace of a telluric star. The effective integration 
time was thus 80 min for the North image, and 40 min 
for the South image. 


4. RESULTS 

4.1. Velocity dispersion of the BCG 

The BCG velocity dispersion ( a v ) was obtained using 
the method of Kelson et al. (2000; see eq. 9 therein). 
Here, redshift and cr v are jointly determined by min- 
imizing x 2 as the galactic spectrum is compared to 
a shifted, broadened, continuum-matched stellar tem- 
plate in pixel-space using a gradient search algorithm 
(e.g., IDL MPFIT). As the SNR of the blue channel BCG 
data was low, we excluded it from this process, using 

12 http: //www. ucolick.org/~xavier/IDL/index.html 


it only to provide initial redshift information based on 
Ca H & K. Using both the blue and red channel spec- 
tra in their entirety, we find a redshift for the BCG of 
zbcg = 0.2373 ± 0.0002. Unfortunately, the loss of the 
blue channel data combined with the BCG redshift and 
the fixed disperser available to remote observers on DIS 
removed any overlap between BCG and observed stellar 
standards, rendering the latter unusable as templates. 
We employed instead a set of high-resolution spectra 


from the UVES Paranal Observatory Project (Bagnulo 


et al.|2003), selected by spectral type without regard tor 


metaihcity. 

Before comparison with the BCG spectrum, it was nec- 
essary to match the template resolution with that of DIS. 
The latter quantity was determined in two ways: 1) gaus- 
sian fitting of sky emission lines on 2-D DIS spectra; 2) 
gaussian broadening of UVES templates and compari- 
son to DIS stellar spectra. The methods yielded values 
of FWHM insi = 1.46 ± 0.15 A and 1.78 ± 0.21 A, re- 
spectively, showing no significant trend with wavelength. 
We adopted the average of these, FWHM inst = 1.62 A 
as the instrumental resolution for DIS, corresponding to 
Ginst ~ 30 km s -1 at the relevant wavelengths, observed 
frame. All templates were broadened to this resolution 
before beginning the BCG fitting procedure. 

At the SNR of our BCG spectrum, the uncertainties as- 
sociated with template variation/mismatching (system- 
atic) and photon statistics (random) are dominant and 
of comparable order. Precisely quantifying the system- 
atic effects requires extensive simulation and is neces- 
sary when investigating detailed physical relations based 
on velocity dispersions, such as Mf, — a . Our interest 
in this quantity, however, extends only to providing a 
prior for the lensing model below. Hence, we deem it 
sufficient to verify that our procedure adequately repro- 
duces accepted values of a v from the literature. We took 
as a test sam ple a se t of spectra for 20 Virgo ellipticals 
from Dressier (1984) for which cr„ measurements were 
made using a region surroun ding the A5176 M g b triplet 
- which we also capture. The Dressier ( 1984| data are of 
comparable SNR to our own. Restricting the fitting to 
the 4900 - 5400 A interval used in the original paper, the 
template that best-reduced \ 2 across all the (Dressier 


1984) galaxies yielded —0.02 < logia(a v /cr V: D 84 ) < 0.05 
for a v> D 84 > 85 km s _1 with a difference in the mean con- 
sistent with zero. Template-to-template variations intro- 
duce a logarithmic scatter of less than 0.1 dex across all 
a v , but do not bias the results. 

Adopting as much of the Dressier (1984) interval as 
our data allow (A G [5024, 5400] A) and masking tel- 
luric features, we find cr„ i bcg = 298 ± 8 (random) ± 
18 (systematic) km s -1 . Values are taken from the best- 
fitting template with the systematic error reflecting RMS 
variation across templates. 


4.2. Lensing model and magnification 

We construct a lens model of the cluster using the 
publicly-available software LENST00L (Julio et al. 2007). 
The cluster and the BCG are each modeled by a 
pseudo-isothermal ellipsoid mass distribution (PIEMD; 
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Limousin et al. 2005rH parameterized by its position 
x , y\ a fiducial velocityaispersion ctpiemd! a core radius 
fcore'i a cut radius r cu t\ ellipticity e = (a 2 — 6 2 )/(a 2 + 6 2 ), 
where a and b are the semi-major and semi-minor axes, 
respectively, and a position angle 9. We fix the BCG 
parameters x, y, e, and 9 at their observed values as de- 
rived from the GALFIT modeling described above, and 
fix the cluster cut radius at 1000 kpc as this cannot be 
constrained by the data. This does not significantly af- 
fect the lensing model. The BCG veloc ity dispersion of 
a v = 298 ± 26 km s -1 measured in i 4.1 is set as a Gaus- 
sian prior. The velocity dispersion of the cluster was not 
measured, but we can estimate it. SDSSJ 1438+1454 
appears in the GMBCG catalog (Hao et al. 2010) with 
a weighted richness Nj^f hted = 9.734; following the cal- 
ibration of Becker et al. (2007), we estimate a velocity 
dispersion of a = 318+ 111 km s -1 . The low richness 
and velocity dispersion indicate that SDSSJ 1438+1454 
is a poor cluster or group, in agreement with the general 
appearance of the cluster in the optical imaging, and the 
small separation between the lensed images of the back- 
ground galaxy. We use this velocity dispersion as a Gaus- 
sian prior. Note that the above velocity dis persio ns are 
converted t o ctpiemd using the relations in (Elfasdottir 
et al. 2007). We also note that the BCG in this cluster 
is very dominant, based on the magnitude difference of 
1.64 mag in the i-band between the BCG and next two 
brightest cluster member s (for typical poor clusters, t his 


difference is < 1.0 mag; [Niederste-Ostholt et al. 2010). 
We therefore limit the cluster center to be not farther 
than 39" from the BCG, which is the separation between 
the BCG and the next brightest galaxy. Dominant sys- 
tems tend to show good alignment between the BCG and 
the cluster halo (Niederste-Ostholt et al. 2010), similar 
to the results from the lens modeling shown below. The 
best-fit model is found via Monte Carlo Markov Chain 
(MCMC) optimization in the image plane, using one to 
four coordinates in each instance of the lensed galaxy as 
constraints (see Figure 6). The best-fit model has an 
image-plane RMS of O'/I. Table 2 lists the best-fit pa- 
rameters and their ler uncertainties. 

Since the lensed images are extended, the magnifica- 
tion is not fixed across different portions of each lensed 
image of the source. The average magnification of each 
lensed image is thus measured as the ratio of areas occu- 
pied by the galaxy in the image plane and source plane. 
In Table 3 we report the average magnification for each 
image of the source. We also report the magnification 
of the core regions of the North and South images, as 
these regions were targeted for spectroscopy (see Figure 
5). Table 3 also includes the total magnification of the 
aggregate of all images of the source; note that the to- 
tal magnification of the source is not a linear sum of the 
magnifications of its instances because the source galaxy 
crosses the source-plane caustics and some of its lensed 
instances (the east and west images) represent only parts 
of the source galaxy rather than all of it. In Figure 7 we 
show the magnification contours of the best-fit model 
and a source-plane reconstruction of the galaxy. To es- 
timate the uncertainty on each magnification, we draw 


13 This profile is formally the same as a dual Pseudo Isothermal 
Elliptical Mass Distribution (dPIE, see Elfasdottir et al. 2007). 


parameters from chains in the MCMC process, compute 
a lens model for each set of parameters, and measure the 
magnification of each image. The uncertainties quoted 
in Table 3 are given as the range of magnifications cor- 
responding to ler in parameter space. 

4.3. Redshift, line fluxes, and metallicity 

We separately fit the Keck NIRSPEC spectra for each 
of the two images of the lensed galaxy, as follows. 
The redshift was set initially by fitting a single Gaus- 
sian to H-a. We then simultaneously fit the [N II] 
doublet and H-a using three Gaussians using the IDL 
Leven berg-Marquardt le ast-squares fitting code MPFIT- 
FUN (Markwardt 2009). Wavelengths were set at the 
NIST values^j and allowed to vary within ±2cr of the 
initial fit to H-a. A common linewidth for all lines was 
allowed to vary freely. The [N II] do ublet flux ratio wa s 
fixed at the value computed in Storey & Zeippen (2000). 

Table 4 reports the redshift for SGAS 1438 as measured 
from simultaneously fitting the [N II] and H-a lines in 
each near-IR spectrum. The result is consistent with 
that found by fitting H-a alone. The weighted mean 
redshift for the two images is z = 0.81589 + 0.00005. The 
two images have identical redshifts within uncertainties. 
The spectra are plotted, with their multicomponent fits, 
in Figure 8. Table 4 lists measured fluxes for H-a and 
[N II], 

We infer metallicities from the measured H-a to [N II 


line flu x ratios using the calibration of |Pettini &; Page 
(2004|. The metallicity is calculated for each image o: 


the lensed galaxy, and tabulated in Table 4. The summed 
[N II] and H-a fluxes over both images yields the most 
precise metallicity: 12+log (O/H) = 8.62+0.06 using the 
linear Pettini & Pagel calibration, and 12 + log (O/H) = 
8.64 ± 0.09 using their third-order polynomial calibra- 
tion. For r eference, sol a r met allicity is 12 + log (O/H) = 
8.69 + 0.05 ( Asplund et al.|2009 ). Thus, we conclude that 
SGAS 1438 has solar metallicity. 

4.4. Spectral Energy Distribution Fits 

We fit the rest-frame UV, optical, and near-infrared 
SED of SGAS 1438 following the procedure outlined in 
( Wuyts e t al.||2010l|2012D. We use Hyperz (|Bolzonella et 
aT 2000) to perform SED fitting at fixed spectroscopic 
redshift with the newest Bruzual & Chariot population 
synthesis models (CB07 - see Bruzual & Chario t (2003) ), 
using the Chabrier initial mass function ( Cha brier] 2003 ) 
and the Calzetti dust extinction law (Calzetti et al. 
2000 ) . We allow a range of exponentially declining star 
formation histories with e-folding times r ranging from 
10 Myr to 2 Gyr as well as a constant star formation 
history. The age of the stellar population is allowed to 
vary between 50 Myr and the age of the Universe at 
z = 0.816(^j Informed by the near-IR spectroscopy, we 
fix the metallicity at the solar value. Figure 9 plots the 
rest-frame UV to near-IR photometry reported in Table 
1 as well as the best-fit stellar population model. From 
this SED, the derived stellar mass, corrected for lensing, 
is log(M*/M 0 ) = 10.8l°J. 

The Herschel photometry shows that SGAS 1438 is 
extremely bright in the far-IR. Its infrared spectral en- 

14 http:/ /www/pa. uky.edu/~peter /atomic 

15 6.556 Gyr in our assumed cosmology. 
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ergy distribution peaks in log (/„) between 90 and 140 
microns in the rest frame, which is a considerably red- 
der peak wavelength than nearby template galaxies like 
M 82 or Arp 220. We fit the 12-5 00 nm phot ometry 
with the local galaxy template set of Rieke et al. ( 2009), 
by convolving each template with the WISE. PACS and 
SPIRE filter curves. As shown in Figure 10, the best- 
fitting template is the log (L(TIR)) = 10.75 L 0 Rieke 
template. (L(TIR) is defined as the rest-frame 8-100 fi m 
luminosity.) Correcting for lensing magnification and in- 
tegrating this best-fitting template from rest-frame 8- 
1000 fim yields an intrinsic far-infrared luminosity of 
L(TIR)= 4.0l^ x 10 11 L 0 . 

It is important to note that the best- fitting template 
was made from the photometry of galaxies almost a 
decade less infrared-luminous than is SGAS 1438: the 
template’s luminosity is 7x lower than our measured 
L(TIR). In other words, templates of comparable lumi- 
nosity provide poor fits, whereas a less luminous tem- 
plate, scaled up in luminosity, provides a much better fit 
to the shape of the far-IR SED. We will return to this 
point in )J5] 

4.5. Star Formation or AGN? 

An important question when deriving physical condi- 
tions for SGAS 1438 is whether there is any significant 
AGN contribution to the observed emission, either from 
an AGN in the lensed galaxy itself, or from an AGN in 
the central BCG (the lens galaxy). Here we review the 
(lack of) evidence for an AGN. The bulk of the light at 
near-IR wavelengths and blueward clearly comes from 
an extended source, both for the lensed source itself and 
the BCG. The DIS optical spectrum of the BCG does 
not show any emission lines: the observed spectra in- 
cludes the expected wavelength of [Oil] 3727. In the 
NIR spectrum of the lensed source, the H-a line is not 
unusually broad or strong. Thus, there is no evidence for 
a significant AGN component in the rest-frame optical or 
near-IR. 

The reasonable match between the FIR SED and the 


star- forming template set of Rieke et al. ( 2009 ) also gives 
some encouragement that any putative AGN contribu- 
tion is no larger than might be disguised in that tem- 
plate set. The far-IR SED is matched by a relatively 
cold template; an AGN would be expected to make hot- 
ter dust and thus shift the peak to the blue, which is 
not obser ved . We note that the Herschel photometry at 
70 fin i (]3.3) also provides a limit on any AGN com- 
ponent from the BCG. The measured flux in the South 
image of SGAS 1438 is 22 ±5% of the total flux. The pre- 
dicted fraction from the lens model is 31^2%, an d from 
the measured flux ratios at shorter wavelengths that re- 
solve all components of the source it is 28 ±2%. Taken at 
face value, the 70 fi m and shorter wavelength flux ratios 
limit the contribution of the BCG at 70 fin i to less than 
21% of the total flux, and are consistent with no contribu- 
tion at all. Finally, there is also no evidence for emission 
at the position of SGAS 1438 in either the FIRST ra- 
dio data (|Becker et al.||1995| or at X-ray wavelengths in 
the ROSA'T All Sky Survey, although these data do not 
set any strong limits on a low-luminosity AGN contri- 
bution. Regardless, absent data at any wavelength that 
suggest the presence of an AGN, we proceed under the 


assumption the emission observed from the lensed source 
is powered by stars only. 

4.6. Star formation rate 

We estimate the star formation rate (SFR) of 
SGAS 1438 from two methods: H-a and the far-infrared. 

To convert the H-a fluxes reported in Table 4 into a 
star formation rate (SFR) captured by t he NI RSpec slit, 
we a pply the calibra tion of |Kennicutt ( 1998) , modified 
for a Chabrier (2003) initial mass function: 


SFR = 0.66 L(Ha) x 7.9 x 10 


-42 


(1) 


With L(Ha ) specified in ergs s -1 this is the apparent 
star formation rate in M 0 yr^ 1 in the slit; it requires an 
aperture correction, magnification, and extinction cor- 
rection to get the intrinsic star formation rate. Poor 
and variable seeing during the Keck observation means 
that slit losses are high and uncertain. We estimate the 
seeing by convolving the J-band image of the brightest 
cluster galaxy until it matches the width of the BCG in 
the two-D spectrum, and then compute the slit losses by 
examining the seeing-degraded images of the source. We 
conclude that the aperture correction is large: a factor of 
4± 1. Final intrinsic SFRs from the H-a measurements 
are computed for each image combining uncertainties in 
the line flux (Table 4), the magnification (Table 3), and 
the aperture correction. We find SFR = 3ll^ M 0 yr 1 
from the North image and SFR = 25+g 6 M 0 yr -1 from 
the South image. Note these are both measurements of 
the total SFR for the entire source; combining with equal 
weights yields a final best value of SFR#_ Q = 28+7 3 
M 0 yr -1 . 

We also estimate the star formation rate from the far- 
infrared luminosity, using the template fitting discussed 
in §4.4| and plotted in Figure 10. The integrated rest- 
frame 8 to 1000 fi m luminosity, corrected by the total 
magnification, is 4.0+}'® x 10 11 L 0 . Via the calibration 
of Kennicutt ( 1998 1 this corresponds to a SFR of 41+}} 
M 0 yr- 1 . 

4.7. Extinction 

If we attribute to extinction the higher value of the 
SFR measured from the FIR versus H-a, then the extinc- 
tion at H-a is 0.4++4 mag. The extinction at H-a esti- 
mated from the SED fit is 0.9+}}} mag. Thus, SGAS 1438 
may have slightly hig her ex tinction than the calibrating 
sample of Kennicutt (1998), though only at ler signifi- 


cance. 

4.8. Physical Size and Structural Properties 

The good image quality of the shorter wavelength 
imaging of SGAS 1438 allows a robust measure of the 
galaxy’s size and morphology. The strong lensing model 
for SGAS 1438 also shows unusually low distortion for 
the images of this source, compared, for example, to the 
highly magnified gian t arcs m o re ty p ical in strong lensing 
galaxy clusters (e.g., Sharon et al. 2012]) ; this suggests 
that measurements in the image plane can be simply 
converted to the source plane. To thoroughly explore 
the source morphology and size, we have fit the North 
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and South images of the galaxy both directly in the im- 
age plane and also in reconstructed images in the source 
plane. We consider the JHK S images as well as the 
i-band image; these are the deepest and best seeing im- 
ages, spanning two distinct instruments and telescopes. 
Visual examination of the initial images of SGAS 1438 
(see Figure 3) suggests a two-arm spiral morphology for 
the source. In fitting the source using simple models we 
mask the obvious portions of these apparent spiral arms. 

To fit in the image plane we use GALFIT and the 
PSF references already established for modeling the BCG 
light, and fit the isolated (and complete) North and 
South images of SGAS 1438 independently in the BCG- 
subtracted frames. Fits were made using both a single 
Sersic profile^] as well as a single Gaussian profile to fa- 
cilitate com parison to the measurements of Rujopakarn| 


formation - the uniform red light seen across much of the 
source is that of a vigorously star forming disk. A possi- 
ble alternate interpretation of the source is as a merger 
with both arms representing tidal tails; blue tidal fea- 
tures with many young star s are seen in mergers nearby 
(e.g., Whitmore et al. 2010). However, the good fit to a 
exponential disk model, and the otherwise uniform col- 
ors of this source disfavor the merger hypothesis. Higher 
resolution imagery will be required to further illuminate 
this question. 

5. DISCUSSION AND CONCLUSIONS 

The far-IR luminosity of SGAS 1438 makes it a “Lu- 
minous Infrared Galaxy” or LIRG. It is instructive to 
place SGAS 1438 in the context of the deep Spitzer 
surveys done using the MIPS 24 /jm band. The ma- 


plane by dividing by the square root of the magnification; 

shifts of 0.6 < z < 1.2 

Le Floc’h et al. 

2005 

). At 

in the case of equal magnification radially and tangen- 

z=U.8 the typical infrared luminosity tor JV 

UPS detec- 


tially, as is approximately true here, this is appropriate. 
The results of this fitting are reported in Table 5. Treat- 
ing each image in each filter with equal weight, the mea- 
sured Sersic index in the image plane is 0.87±0.14 with 
an effective radius of 4.1±0.6 kpc and a Gaussian size 
(in cr) of 7.7±0.8 kpc. To measure morphology in the 
source plane we construct source plane images from the 
image plane using the lensing model, by simple interpola- 
tion on the irregular grid of image pixels produced in the 
source plane. This is not optimal from a signal to noise 
perspective since it treats all image plane pixels equally, 
regardless of magnification; it is however sufficient for 
the purposes of this paper. To obtain a PSF reference, 
we map the image plane PSFs back to the source plane 
in the same way, at the position of the center of the 
source. Formally, in the source, the realized PSF is a 
function of position; again, because the distortions are 
not strong, this effect is sub-dominant. We have tested a 
range of PSFs sampled from locations around the center 
of the source, and find no significant effect on the re- 
sults. Again treating each image in each filter with equal 
weight, the measured Sersic index in the source plane is 
1.23±0.25 with an effective radius of 4.3±0.6 kpc and a 
Gaussian size of 7.1L0.9 kpc. 

These measurements, which are consistent across four 
filters, two images of the source, and both the image 
plane and the source plane, clearly show that SGAS 1438 
is dominated by a disk with an exponential light distribu- 
tion. This is consistent with the visual impression of the 
system as a two arm “grand design” spiral galaxy. The 
nucleus of SGAS 1438 also appears significantly elon- 
gated, suggesting a bar. 

4.9. Color and color gradients 

Figures 6 and 7 include both color composite images as 
well as greyscale images reconstructed in the source plane 
in various filters. The images show that SGAS 1438 has 
an overall red and uniform morphology. The most sig- 
nificant deviations from color uniformity occurs in one of 
the putative spiral arms, which shows a few bluer clumps. 
Note that only the core of SGAS 1438was targeted for 
near-IR spectroscopy and shown to have vigorous star 

16 An exponential disk has a Sersic index of 1, and a de Vau- 
couleur profile has a Sersic index of 4. 


tions is L(TIR)^ 1 x 1 0 11 Lfp, with uncertai nties of or- 


der of a factor of 2-3 (Le Floc’h et al. 
shift and luminosity of SGAS 1438 o 


2005). The red- 
z = 0.816 and 


L(TIR)= 4.0 x 10 11 Lq mean that were SGAS 1438 not 
lensed, it would appear as a rather typical actively star- 
forming galaxy in the Spitzer deep fields. Such galaxies 
dominate the star formation rate density of the universe 
a t that epoch (|Pa povich et al. 2004 Le Floc’h et al. 2005 1 


Rujopakarn et al. 1 2012 1 succinctly summarize results 
from the past decade showing that actively star-forming 
galaxies in the distant universe are different from z = 0 
analogues of comparable luminosity in at least three ma- 
jor ways. First, their infrared SEDs are colder. Sec- 
ond, their aromatic features grow stronger with redshift 
for fixed L(TIR). Third, their star formation occurs on 
physically larger scales. 

The first two effects were discovered first, and a com- 
mon cause was suspected: that if star formation was 
happening on large spatial scales, that would naturally 
explain the higher aromatic feature strengths (less ex- 
tinction) and colder SEDs (less intense radiation fields). 
New mea surements of large sizes for z > 1 LIRGs and 
ULIRGs (Rujopakarn et al. 20111 confirm those predic- 
tions. Thus, we have a framework to understand how 
the population of vigorously star-forming galaxies evolve 
with redshift. Though SGAS 1438 is only one galaxy, its 
high lensing magnification enables us to test this frame- 
work at otherwise unavailable photometric quality and 
spatial resolution. Thus, we now summarize the evidence 
for evolution with redshift in the SED shape, aromatic 
strength, and size of star-forming regions in luminous in- 
frared galaxies, and discuss how SGAS 1438 tests this 
picture. 

Cooler SEDs: That luminous infrared galaxies in the 
distant univers e may have co l der SED s t han a t z = 0 was 
first noted by Rowan-Robinson et al. (2004) from ISO 


data for LIRGs at 0.15 < z < O.J^ and then confirmed 


with Spitzer (Rowan-Robinso n et al. 2005 Sajina et al. 
2006| |Symeonidis et ai.||2009p. Spitzer 70 and 160 //m 
photometry showed that LIRGs and ULIRGs at z ~ 1 
have, on average, colder far-infrared SEDs th at peak at 
longer wav elengths than z ~ 0 analogues ( Syme onidis 
et al.||2009|). Figure 5 of |Symeonidis et al.|(2009|) illus- 
trates how 0.1 < z < 1 LIRGs have longer-wavelength 
SED peaks (in /„) than z = 0 LIRGs. The SED peak of 
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SGAS 1438 is entirely consistent with their distant LIRG 
sample, and inconsistent with the z = 0 sample. Those 
authors note that a simple way to describe the evolution 
in the luminosity-SED relation with redshift is that “the 
far-IR SEDs of many high-redshift sources may resemble 
more closely those of lower-luminosity galaxies locally.” 
This result was better qu antified wit h m ultiband pho- 


tomet ry fro m Herschel: Elbaz et al. (2010) and Hwang 
et al. (2010) found that the median dust temperature of 
LIRGs at 2 ~ 1 is several K colder than z = 0 analogues; 
the offset is even more pronounced at higher luminosities. 

The dust-reprocessed SED of SGAS 1438 is entirely 
consistent with these trends. We fit the SED with the 


local template set of Rieke et al. (2009), with the lu- 
minosity scaling as a free parameter. The best-fitting 
template belongs to z = 0 galaxies that are nearly an or- 
der of magnitude less IR-luminous than SGAS 1438. The 
dust SED of SGAS 1438 looks like a cool, low-luminosity 
galaxy that has been “scaled up” in luminosity. 

Aromatic feature strength: The evolution toward 
stronger aromatic features with redshift, at fixed L(TIR), 
was first discovered from Spitzer mid-infrared spectra 
and far-I R spectral energy distributions of 1 < z < 3 


galaxies dPapovich et al. 2007 Rigby et al . 2008 Far- 
rah et al.j 2008 Mencndez-Dcimestre et al. 2009 ) . The 
explanation proposed by the latter three sets of authors 
was that star formation in distant galaxies may occur on 
larger spatial scales than in local galaxies. Unfortunately, 
for SGAS 1438 we cannot isolate the aromatic feature 
strengths, as we lack mid-infrared spectra for this source. 
The WISE photometry at 12 and 22 /rm do contain flux 
from the aromatic features, but blended in unknown ra- 
tio with hot dust continuum. Thus, it is not possible to 
test whether the aromatic feature s in SG AS 14 38 ob ey 
the evolutionary relation noted in iRigby et al. (2008). 


Size: Rujopakarn et al. (201 If measured Gaussian 


sizes of 0.4 < z < 2.5 LIRGs and ULIRGs of 2-10 kpc. 
These are considerably larger than the the measured 
< 1 kpc sizes of star form ation in comparable-lumin osity 
galaxies at z = 0. Thus, Rujopakarn et al. (2011) con- 
clude that luminous star formation episodes at z = 0 are 
found only on sub-galactic (D < 1 kpc) spatial scales, 
while at z > 0.4 they occur on larger, galaxy-wide scales 
( D ~ 3-8 kpc), with consequently lower star formation 
rate surface densities. 

The Gaussian size of SGAS 1438 (7.4±0.6 kpc) is 
slightly larger than the sizes of 3-6 kpc measured by 
Rujopakarn et al. (2011) for z ~ 1 LIRGs, but consis- 


tent given our uncertainties and their small sample size. 
We can also compare to size measurements from HS T 
ima ging o f disk galaxies in the GEMS fields. Barden et 


al. (2005) predict a relationship between effective radius 
{of a Sersic profile fit) and total stellar mass. For the 
measured stellar mass of SGAS 1438 of ~ 6 x 10 10 M 0 , 
the predicted effective radius from their mean relation 
is ~ 7 kpc. This is somewhat larger than the measured 
value of 4.3±0.4 kpc for SGAS 1438, though consistent 
with expectations given the significant spread in the val- 
ues of effective ra dius at fixed st ellar mass apparent in 


the results of [Barden et al. (]2005). 

SGAS 1438 in this context: SGAS 1438 is only 
a single galaxy, but the strong lensing magnification al- 
lows for a detailed suite of measurements otherwise im- 
possible to obtain. The redshift, luminosity, and far-IR 
SED of SGAS 1438 show it to be a typical vigorously 
star forming galaxy at z ~ 0.8, and thus representative 
of those galaxies that contribute most of the universe’s 
star formation at this earlier epoch. Its dust-reprocessed 
SED is clearly cooler than comparably luminous galaxies 
at z = 0, and comparable to its peers at similar red- 
shift. The spatial extent of star formation in SGAS 1438 
is much larger than in comparably luminous galaxies at 
z = 0, and comparable to (or slightly larger than) the 
half-dozen LIRGs measured at 0.5 < 2 < 1 by |Ru-| 


jopakarn et al. (2011), confirming that result at higher 
signal-to-noise and spatial resolution. The observations 
of SGAS 1438 presented here knit together both the size 
and SED-to-luminosity evolution seen piecemeal in other 
work in a single, and apparently typical, object. 

What is unique about SGAS 1438 is that lensing mag- 
nification allows us to determine the morphology of this 
vigorous star formation at high spatial resolution and 
signal-to-noise ratio. Morphological measurements, both 
in the image plane and the reconstructed source plane, 
strongly suggest that the star formation in SGAS 1438 
is occurring in a large exponential disk. The surpris- 
ing result, as show in Figure 7, is that star formation 
appears to be distributed reasonably smoothly over the 
entire disk of this galaxy. 

Thus, we confirm previous results that the mode of 
star formation is different at z ~ 1 than at z = 0. Fur- 
ther, at high spatial resolution we image that different 
mode, and see that in SGAS 1438, vigorous star for- 
mation is taking place smoothly across an exponential 
disk with an effective radius of 4.3±0.6 kpc and an ap- 
pearance suggestive of a two-armed ’grand design’ spiral 
morphology. 


APPENDIX 
THE RED HERRING 

The discovery of cluster SDSSJ 1438+1454 as a strong lens has a curious history, in that the initially-identified arc 
candidate proved not to be lensed at all, while the follow-up infrared imagery revealed a different, very red galaxy 
that was also suggested in the SDSS imaging, is actually lensed, and is the subject of this paper. For completeness, 
we now describe the initial arc candidate. 

The lensing search process detailed in Gladders et al. (in prep.) does not record which feature caused a particular 
field to be flagged as a candidate for lensing; the most likely feature is this arclet-like source, though a second faint and 
red feature - the lensed source discussed in this paper - is just visible in the SDSS imaging. A portion of the custom 
SDSS color image used in the lensing search of Gladders et al. is reproduced in Figure 1. 

Follow-up imaging from the Nordic Optical Telescope (Figure 1) strengthened the lensing interpretation for the 
candidate arclet, since the curvature suggested in the SDSS imaging is confirmed, and the image does not show any 
bulge-like component which might be expected in the case of an edge-on galaxy. 
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The candidate blue arclet was also observed spectroscopically using the ARC 3.5m telescope and the Double Imaging 
Spectrograph (DIS) on March 2nd, 2008. These data show that the candidate blue arclet is not in fact lensed; rather, 
it has a redshift of z = 0.231 from clear detections of [Oil] 3727 and H -a, which is consistent with this source being 
a cluster member, given the SDSS redshift for a nearby apparent cluster galaxy at z = 0.237 as well as the redshift of 
the BCG of 0.2373 (see §4.1). 
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TABLE 1 

Imaging observations and photometry of SGAS 1438 


Telescope &; Instrument 

filter 

int (s) 

PSF (") 

date 

m ab 

Gemini-N GMOS 

SDSS g 

300 

0.61 

29 Apr 2011 

21.45±0.16 

Gemini-N GMOS 

SDSS r 

300 

0.56 

29 Apr 2011 

20.24i0.il 

Gemini-N GMOS 

SDSS i 

300 

0.54 

29 Apr 2011 

18.85±0.10 

Gemini-N GMOS 

SDSS z 

300 

0.50 

29 Apr 2011 

18.10±0.09 

Magellan FourStar 

Y 

800 

0.68 

19 Feb 2011 

18.03±0.12 

Magellan FourStar 

J 

800 

0.66 

19 Feb 2011 

17.34±0.11 

Magellan FourStar 

H 

800 

0.63 

19 Feb 2011 

16.84±0.11 

Magellan FourStar 

K s 

1572 

0.53 

23 Jan 2011 

16.31±0.11 

Spitzer IRAC a 

3.6 pm 

600 

1.7 

19 Sep 2010 

15.74±0.13 

Spitzer IRAC a 

4.5 pm 

600 

1.7 

19 Sep 2010 

16.14±0.13 

WISE 

12 pn i 

167 

6.5 

Jan/Jul 2010 

14.96±0.08 c 

WISE 

22 pm 

193 

12.0 

Jan/Jul 2010 

14.10±0.21 c 

Herschel PACS b 

70 pm 

900 

5.6 

09 Jul 2011 

12.06±0.15 

Herschel PACS b 

100 pm 

900 

6.8 

09 Jul 2011 

11.11±0.15 

Herschel PACS b 

160 pm 

1800 

11 

09 Jul 2011 

10.00±0.15 

Herschel SPIRE b 

250 pm 

185 

18 

12 Jul 2011 

10.06±0.15 

Herschel SPIRE b 

350 pm 

185 

25 

12 Jul 2011 

10.41±0.15 

Herschel SPIRE b 

500 pm 

185 

36 

12 Jul 2011 

11.34±0.15 


Note. — Columns are: 1) Telescope and Instrument; 2) filter or bandpass; 3) integration 
time in seconds; 4) point spread function, quoted as the seeing disk full width at half maximum 
(FWHM) for ground-based instruments, and the tabulated point spread function FWHM for 
space observatories; 5) universal date of observation; 6) AB magnitude for all images of the 
source in aggregate - see main text for further details; 

a Spitzer data are from Program 70154, “Mass across the Redshift Desert: Stellar Masses in a 
Large and Uniform Sample of Strongly-Lensed Galaxies at 1 < z < 3” (PI Gladders). 
k Herschel images are from program DDT_jrigby_2, “The far-IR SED of a highly magnified red 
galaxy” (PI Rigby). 

c Central galaxy flux subtracted using SED model, as detailed in main text. 


TABLE 2 

Best-fit lens model parameters 


Halo 

RA 

Dec 

e 

e 

Tcore 

f'cut 

& PIEMD 

(PIEMD) 

(") 

(") 


(deg) 

(kpc) 

(kpc) 

(km s -1 ) 

Cluster 

ij 7+2.3 
' -'-4.9 

9 9+1-6 
z>z -0.8 

to o 

O CO 

o d 
+ 1 

LO 

d> 

165^ 


[1000] 

672^222 

BCG 

[0] 

[0] 

[0.64] 

[167] 

i o+2. 0 
L ‘ 6 -1.0 

cy+23 
D ' -46 

267 +71 

-30 


Note. — All coordinates are measured in arcseconds relative to the center of 
the BCG, at [RA, Dec] = [219.68768 14.903507]. The ellipticity is expressed as e = 
(a 2 — 6 2 )/(a 2 + 6 2 ). 6 is measured north of West. Error bars correspond to 1 cr 
confidence level as inferred from the MCMC optimization. Values in square brackets 
are for parameters that are not optimized. 
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TABLE 3 

Magnifications of images of 
SGAS 143845.1+145407 


Image 

Magnification 

Notes 

North 

6-5+j 

complete image 

South 

5-2tk? 

complete image 

East 

z ' y -0.8 

partial image 

West 

4 n+ 2 - 5 
^• u — 1.3 

partial image 

North-Core 

4-6+1 

NIR slit region 

South-Core 

5-8+2 

NIR slit region 

All 

16-8+2 



Note. — Columns are: 1) a source image 
name - see Figure 6 2) the estimated magnifica- 
tion from the strong lensing model, and 3) notes 
about the nature of each image. Due to the lens- 
ing configuration two images are complete im- 
ages of the source, and the other two are partial 
images only. The magnification of the ‘All’ im- 
age gives the magnification of all four images in 
aggregate relative to a single unlensed instance 
of the source; this is the appropriate magnifica- 
tion to use when converting the observed pho- 
tometry in Table 1 back to unlensed quantities. 


TABLE 4 

Measured emission lines in SGAS 143845.1+145407 


object 

redshift 

f(H-o) 

f(NII 6585) 

[N IIJ/H-q 

12±log(0/H) 

12±log(0/H) 

North image 
South image 
both images 

0.81588 ± 0.00006 
0.8159 ± 0.0001 
0.81589 ± 0.00005 

229 ± 16 
238 ± 28 
467 ± 32 

83 ± 14 
69 ± 20 
152 ± 24 

0.36 ±0.07 
0.29 ±0.12 
0.33 ±0.07 

o fic+0.04 

0 -O 0 _o 05 

o cq+ 0.07 
8.oy_0 09 
o f?o+0.05 
0.06 

° ,by -0.09 

8 59 +0 ' 12 

°- oy -0.14 
o £/j+0.09 
°- b4 -0.09 


Note. — Columns are: 1) object name, 2) measured redshift and uncertainty from the three-component 
fit to H-a and [N II], 3) H-o: flux and uncertainty, 4) [N II] 6585 A flux and uncertainty, 5) ratio of [N 
II] /H a and uncertainty, 6) derived oxygen abundance using the linear N2/H-a calibration of Pettini & 
Pagel 2004; 7) derived oxygen abundance using the third-order N2/H-a calibration of Pettini & Pagel 
2004. Reported fluxes are in units of 10 — 17 erg s -1 cm -2 . We separately quote measured values for the 
North and South images. Where we quote values for “both images”, these are the weighted mean and 
error in the mean over all components for the redshift, the sum over all components for each flux, and the 
ratio of the summed fluxes for the line ratio. 


TABLE 5 

Morphology and Size of SGAS 143845.1+145407, Sersic and Gaussian Profile Fits 


Image 

Filter 

R e (kpc) 

Image Plane 

Sersic Index Gaussian a (kpc) 

Re (kpc) 

Source Plane 

Sersic Index Gaussian cr (kpc) 

South Image 

K s 

3.6 

0.99 

7.0 

3.7 

1.19 

7.0 


H 

3.6 

0.96 

7.2 

4.6 

1.50 

8.2 


J 

4.0 

0.71 

8.4 

4.8 

1.66 

8.4 


SDSS-i 

4.3 

0.74 

7.7 

5.1 

1.22 

7.0 

North Image 

K s 

4.0 

1.01 

6.9 

3.9 

1.18 

6.1 


H 

4.1 

1.01 

7.3 

3.8 

1.20 

6.2 


J 

4.9 

0.73 

9.1 

3.7 

1.03 

6.4 


SDSS-i 

4.6 

0.78 

7.8 

4.5 

0.87 

7.2 

Average 

all 

4.1±0.6 

0.87±0.14 

7.7±0.8 

4.3±0.6 

1.23±0.25 

7.1±0.9 


Note. - — Columns are: 1) source image name 2) filter, 3)-5) effective radius, Sersic index and Gaussian size in the 
source plane from fits in the image plane, and 6)-8) effective radius, Sersic index and Gaussian size from fits to the 
reconstructed source plane images. In the case of the image plane, effective radii and Gaussian sizes are computed in 
the source plane simply by dividing by the square root of the magnification, as this is the limiting case suggested by 
the only modest distortions seen in the magnified images in the image plane. No such correction is applied to fits in 
the source plane, however the treatment of the PSF is commensurately more complex, as detailed in the main text. 
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Fig. 1. — Left Panel: A 90x90” image cutout from SDSS DR7 imaging data (blue: g-band, green: r-band, red : i-band + z-band). The 
faint red light around the BCG, which is visible all but to the right of the BCG, is the lensed source discussed in this paper. The arclet to 
the lower left is the candidate lensed source discussed in Appendix A. Right Panel: A 45x45 // image cutout from the NOT 2.5m b-band 
image, confirming the curvature of the apparent arclet candidate, as discussed in Appendix A. 



Fig. 2. — A 45x45 // cutout from the Spitzer IRAC 3.6 fim image (greyscale), with the NOT g-band image over-plotted (red contours). 
Two extremely red sources above and below the BCG are immediately apparent; these are two images of SGAS 1438. 
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Fig. 3. — 45x45” cutouts in the grizYJH and K s bands. All images have been scaled to approximately the same greyscale level for the 
BCG. Note the obvious multiplicity of images, particularly in the near-IR, which confirms the lensing interpretation. 


g-band 


i-band 



Ks-band 




Fig. 4. — The images after subtraction of the brightest cluster galaxy (BCG). Residuals near the center of the field are modest, with 
amplitudes that are 1-2% of the BCG core prior to subtraction. These images are rotated by 17 degrees counter-clockwise relative to the 
other on-sky figures in this papaer. 
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Fig. 5. — A 45x45” K s -band image showing how the Keck NIRSPEC slit was placed to cover two images of the lensed source. The slit 
is placed at one of the two nod positions. 



Fig. 6. — The lens model: critical curves (red lines) and caustics (yellow lines) from the best-fit lens model are overplotted on a color 
composite image (red: J, green: i, blue: g ). The coordinates within each arc that were used as lensing constraints are marked in cyan. 
The location of the source in the source plane (see Figure [7]) is indicated by a white box. 
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Fig. 7. — The lens model: left: Contours of absolute magnification are overplotted on a color composite image (red: J, green: i , blue: 
g ) of SGAS 1438. Right: Reconstruction of the source of SGAS 1438, from the best-fit lens model. We ray-trace each pixel in the image 
plane through the lens model, and resample the source plane into a C//05 pixel grid. 



Fig. 8. — Keck/NIRSPEC spectra of Ha and [N II] in SGAS 1438. Spectra of the North and South images are plotted in the left and 
right panels, respectively. Each X-axis shows observed wavelength in fi m, and each Y-axis shows observed specific intensity f\ in units 
of 10 — 17 erg s — 1 cm -2 A -1 . In each panel, the spectrum is plotted in black, the la uncertainty spectrum is plotted in grey, and the 
multi-component Gaussian fit is plotted in blue. The extent of the blue line shows the region over which the continuum was fit. 
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1 5 

Observed wavelength [fim] 


Fig. 9.— The photospheric spectral energy distribution. Datapoints correspond to gri photometry from Gemini, Y JHK S photometry 
from Magellan, and 3. 6-4. 5 //m photometry from Spitzer, as detailed in Table 1. Also plotted is the best-fitting Bruzual Sz Chariot stellar 
population synthesis model, as described in §4.4| The units of the Y axis are AB magnitudes. 
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10 100 
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Fig. 10. — The dust-reprocessed spectral energy distribution. Datapoints correspond to photometry at observed wavelengths of 12 and 
22 fim from WISE, at 70, 100, and 160 fim from Herschel/PACS, and at 250, 350, and 500 ^m from Herschel/SPIRE, as detailed in Table 
1. W aveleng t h ranges fo r each observation are plotted toward the bottom. Also plotted is the best-fitting local— galaxy template (blue line) 
from|Rieke et al.| (j2009b , which is the log (L(TIR)) = 10.75 L© template. We measure the L(TIR) of SGAS 1438 as 7x higher than this 
fiducial template; m other words, the template that best fits SGAS 1438 was made from relatively low-luminosity galaxies, and has been 
scaled up in luminosity by a factor of 7 to fit a galaxy with a much higher luminosity. To illustrate the quality of the template match, we 
also plot the next-higher and next-lower luminosity Rieke templates at log (L(TIR)) = 11.0 L© and log (L(TIR)) = 10.5 L© respectively 
(dashed lines). 


